Quantum measurement corrections to chemically induced dynamic nuclear 

polarization 
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Chemically induced dynamic nuclear polarization has emerged as a universal signature of spin 
order in photosynthetic reaction centers. Such polarization, significantly enhanced above thermal 
equilibrium, is known to result from the nuclear spin sorting inherent in the radical pair mechanism 
underlying long-lived charge-separated states in photosynthetic reaction centers. We will here show 
that the recently understood fundamental quantum dynamics of radical-ion-pair reactions open 
up a new and completely unexpected venue towards obtaining CIDNP signals. The fundamental 
decoherence mechanism inherent in the recombination process of radical pairs is shown to produce 
nuclear spin polarizations on the order of 10 4 times or more higher than thermal equilibrium values at 
low fields relevant to natural photosynthesis in earth's magnetic field. This opens up the possibility 
of a fundamentally new exploration of the biological significance of high nuclear polarizations in 
photosynthesis. 



The quantum dynamics of photosynthesis have, quite 
naturally, attracted a lot of attention recently. Un- 
derstanding if and how Nature exploits quantum 
(de) coherence would have tremendous scientific and tech- 
nological impact. For what follows, photosynthesis can 
be simply visualized as a two-step process: (i) light har- 
vesting and (ii) photochemical processing. Light har- 
vesting molecules absorb the incident photon and guide 
the subsequent excitation to the photosynthetic reaction 
center (RC). It is there where the photochemistry takes 
place, transforming the initial electronic excitation in a 
transmembrane proton pump that further drives the life- 
sustaining chemistry. Quantum coherence effects in light 
harvesting complexes (time scale of 1-500 fs) have been 
addressed theoretically and experimentally [1-9]. We will 
here focus on RC quantum dynamics (time scale of lps - 
10 fis) schematically depicted in Fig. 1A. Photosynthetic 
reaction centers exhibit a cascade of electron transfer re- 
actions, starting from a photo excited donor-acceptor 
dyad PI (chlorophyl-type molecules) and shelving the 
electron further apart to quinones (Qa, Qb), producing a 
long-lived charge-separated state driving the transmem- 
brane proton pump. In each of those steps there is a 
radical-ion pair formed, further perplexing the dynam- 
ics, as radical-ion pairs (RPs) exhibit intricate magnetic 
field effects [l(| . Radical-ion pairs [HI E2 are biomolec- 
ular ions with two unpaired electrons (denoted by the 
two dots in Fig. 1) and any number of magnetic nuclei. 
The magnetic nuclei of the donor and acceptor molecules 
couple to the two electrons via the hyperfine interaction, 
leading to singlet-triplet (S-T) mixing, i.e. a coherent 
oscillation of the total electron spin state, also affected 
by the electrons' Zeeman interaction with the external 
magnetic field. A spin-dependent charge recombination 
terminates the coherent S-T mixing and leads to the reac- 
tion products, either the original neutral dyad PI (singlet 
product) or a triplet intermediate ( 3 PI). Singlet (triplet) 
radical pairs recombine at the singlet (triplet) recombi- 
nation rate, denoted by ks (kx). Triplet products may 
close the reaction cycle by intersystem-crossing into the 
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FIG. 1: (A) Cascade of electron transfer reactions taking 
place in photosynthetic reaction centers. (B) Each charge- 
separated state in (A) is a radical-ion pair which undergoes 
additional quantum dynamics due to internal magnetic inter- 
actions, embodied by the Hamiltonian 7-L, of the two unpaired 
electrons with each-other, the magnetic nuclei and the exter- 
nal magnetic field. These magnetic interactions can give rise 
to singlet-triplet conversion during the finite lifetime of the 
charge-separated state, further perplexing the dynamics of 
the whole reaction. Figure adapted from the review of Boxer, 
Chidsey and Roelofs [ly|. 



It is apparent from Fig. 1 that when the lifetime of 
the charge-separated states (which increases with sepa- 
ration) is such that magnetic-field effects (with typical 
timescale of e.g. S-T mixing on the order of 0.1-1 fis) 
can come into play, radical-pair dynamics get convolved 
into the reaction center dynamics. One aspect of this 
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interrelation pervading a large number of photosynthetic 
reaction centers is chemically induced dynamic nuclear 
polarization [13-20], the enhancement (by several orders 
of magnitude) of the ground state (PI) nuclear spin polar- 
ization resulting from the RP spin dynamics. CIDNP has 
emerged as a rather universal signature of spin order in 
photosynthetic reaction centers [2l| and its possible oper- 
ational significance for photosynthesis is an open funda- 
mental question. CIDNP is based on nuclear spin sorting, 
to be shortly explained, taking place in radical-pair reac- 
tions. Three fundamental mechanisms have so far been 
known to produce CIDNP signals, the three-spin-mixing 
(TSM) the differential decay (DD) and the differ- 
ential relaxation (DR) mechanism (23[. Note that the 
first is a coherent mechanism resulting from particular 
magnetic interactions within the RP, while the other two 
are incoherent mechanisms, of which the second (DR) 
will not be relevant for the following discussion. We will 
here show that there is yet another, completely unantic- 
ipated mechanism at work at the fundamental quantum 
dynamical level of RP reactions. It is due to the inher- 
ent quantum measurement of the RP's spin state con- 
tinuously going on in the RP as part of its very charge 
recombination process. The decoherence caused by this 
" internal" quantum measurement is responsible for a new 
kind of nuclear spin sorting, particularly efficient at low 
fields, and completely unrelated to other low-field CIDNP 
studies [13]. This new mechanism is ubiquitous, as it is 
largely independent of the particular magnetic interac- 
tions in the radical pairs, and thus opens a new venue 
towards exploring the biological significance of CIDNP- 
induced spin order in natural (earth's field) photosynthe- 
sis. In the following we will (i) briefly explain the basics 
of CIDNP signal generation for the general reader and 
(ii) recapitulate the recent progress in understanding the 
fundamental quantum dynamics of radical-ion-pair reac- 
tions. Combining (i) and (ii), we will then explain the 
quantum measurement corrections to CIDNP. 

Consider a radical-ion pair with just one nucleus cou- 
pled to the electron of e.g. the donor. At time t = the 
radical-ion pairs start out in the electronic singlet state, 
while the nuclear spins are thermally polarized, which for 
protons at earth's field corresponds to a polarization of 
about 10~ 10 , practically taken to equal zero. Thus, in 
the ensemble of radical-ion pairs half will have the nu- 
clear spin in the |+) state, while in the other half the 
nuclear spin will initially be in the |— ) state. The lo- 
cal magnetic field seen by the donor electron will be the 
external field plus the nuclear magnetic field, while the 
acceptor electron will only experience the external mag- 
netic field. Consequently, in the first half radical pairs 
the Larmor precession frequency of the donor electron, 
and thus the S-T mixing frequency, will be slightly higher 
than in the other half (due to the opposite contribution 
of the nuclear spin), hence the coherently formed triplet 
radical pairs will have their nuclear spins predominantly 
in the |+) state, leaving the singlet radical pairs with 
an equal and opposite nuclear polarization. While the 



nuclear spins are thus sorted among singlet and triplet 
RPs, and unless there is a coherent mechanism like TSM 
producing net nuclear polarization, the expectation value 
(I z ) = Tr{pl z } of the z-component of the single nuclear 
spin under consideration is still zero. Here p is the spin 
density matrix of the radical pair, describing the electron 
and nuclear spin degrees of freedom altogether. Accord- 
ing to the traditional understanding of the radical-pair 
mechanism and the concomitant CIDNP signal genera- 
tion, for this nuclear spin sorting to result in a diamag- 
netic ground state with net nuclear polarization, and in 
the absence of any coherent mechanism like TSM, at least 
one of the other two aforementioned incoherent mech- 
anisms should be at play in order to observe nonzero 
CIDNP signals. Here we are going to show that this is 
not the case. To this end we will consider a magnetic 
Hamiltonian that does not exhibit the TSM mechanism, 
and we will also take kg — kx, i.e. we consider the singlet 
and triplet recombination rates of the RP to be equal, 
hence the differential decay (DD) mechanism is not op- 
erational. 

To illustrate the above considerations, we consider a 
radical pair with one spin- 1/2 nucleus isotropically cou- 
pled to the donor electron. In Fig. 2A we display the 
undisturbed spin state evolution brought about by the 
magnetic Hamiltonian TL = Al ■ s± + u>(s± z + S2 Z ), where 
uj is the electron spin Larmor precession frequency and 
A the donor electron hyperfine coupling with the donor 
nuclear spin (in units of frequency). What is plotted in 
Fig. 2A is the spin state's singlet character, given by 
the expectation value of the singlet projection operator, 
(Qs) = Tr{pQs}, and assuming an RP with infinite life- 
time. The nuclear spin, J 2 , of singlet and triplet RPs, 
given by (I z ) s = (QsLQs) and (I Z ) T = {Q t I z Qt), is 
shown in Fig. 2B and 2C, respectively. The previous 
considerations regarding nuclear spin sorting are easily 
visualized by noting (dashed line Li in Figs. 2B and 
2C) that although (I z ) s ^ and (I Z ) T ^ 0, the nu- 
clear spin oscillations in the singlet and triplet subspace 
are exactly out of phase, so the net RP nuclear polar- 
ization (I z ) s + (I Z ) T = 0. This is expected since the 
particular Hamiltonian % does not qualify for the TSM 
mechanism [221 ] . while the differential decay mechanism 
(DD) involves the recombination reaction of RPs and so 
far is out of the discussion. Now, besides the coherent 
RP dynamics driven by dp/dt = — i[H,p] and produc- 
ing the oscillations in Figs. 2A-2C, the spin-selective 
charge recombination of RPs has to be accounted for. 
This is done by a stochastic Liouville equation of the 
form dp/dt — —i[H,p] + C(p), where C(p) is the reac- 
tion super-operator, written in terms of the singlet and 
triplet projection operators, Qs and Qt, and the singlet 
and triplet recombination rates, ks and kx, respectively. 
Recently, quantum measurement theory has been shown 
by Kominis [25], [26[ to fundamentally describe the reac- 
tion dynamics of RPs, resulting in a new form of C(p), 
conceptually departing from the traditional theory used 
until now [271 ]. The main physical point of departure 
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FIG. 2: Time evolution of (A) (Qs), (B) the nuclear spin 
of singlet RPs (QsIzQs) and (C) the nuclear spin of triplet 
RPs (QtIzQt) , when the radical-pair spin state undergoes 
unceasing coherent S-T mixing driven by dp/dt = — i[H, p], 
where H = AI -Si + w(si z + S2 Z ). For this example it was u = 
A/10. The dashed line Li is a reminder of the fact that the 
singlet and triplet RP nuclear polarizations are exactly out of 
phase and sum up to zero. (D) Evolution of (I z ) according to 
the full density matrix equation of Kominis [26|], with ks = 
kr ~ 4A. The absolute value of the nuclear polarization is 
shown in the right y-axis, while its enhancement relative to 
the thermal equilibrium polarization, on the order of 10 4 , is 
shown in the left y-axis. Qualitatively the same result (only 
higher by about a factor of 2) is obtained from the Jones- 
Hore theory 30]. In contrast, the traditional theory predicts 
(I z ) = 0, exactly the same as the sum of (B) and (C) that 
applies to the case of no reaction. Note the different time 
scales in (A)-(C) and (D): since k ^> lj, the build-up of (I z )qc 
in (D) takes place at the very beginning of the S-T mixing 
oscillations shown in (A)-(C). 



is the fact that the very charge recombination process 
of RPs damps their coherent S-T mixing, i.e. it brings 
about an unavoidable intra-molecule decoherence in their 
spin state evolution. The traditional theory for the reac- 
tion operator of RPs, although successful in accounting 
for a large volume of experimental data, does not em- 
body this decoherence process, neglecting which leads to 
several conceptual problems [28|, |29[. In any case, in the 
example treated above, if the nuclear spin expectation 
value is calculated according to the full master equation 
of the traditional theory, taking ks = kr, it is found to 
be identically equal to zero, exactly as it is if only the 
Hamiltonian evolution is taken into account. According 
to the common wisdom, this should be perceived as an 
obvious result, since, with ks = &t the differential decay 
mechanism (DD) cannot lead to any net nuclear polar- 
ization. 

However, and this is the counter-intuitive new result 
presented here, if the new master equation derived by 
Kominis (2|| is used, it is found that there actually is a 
non-zero RP nuclear spin polarization depicted in Fig. 
2D. A similar result, different by just a factor of 2, fol- 
lows from another variant of the quantum-measurement- 
based stochastic Liouville equation by Jones and Hore 
[30l | . In this particular example, this nuclear polarization 
is about 10 4 times higher than the thermal equilibrium 
proton polarization (left y-axis of Fig. 2D) and roughly 
independent of ui as long as lj is relatively small (to be 
further discussed later) . What is the source of this new 
kind of CIDNP signal ? We will now go beyond the par- 
ticular example considered above and provide the expla- 
nation for the general case of this new channel for pho- 
tochemically enhancing nuclear polarizations: the funda- 
mental singlet-triplet decoherence underlying the radical- 
pair mechanism. 

To visualize this mechanism, suppose that at time t the 
single-nuclear-spin RP ensemble is in a state described 
by the density matrix p. Assuming a magnetic Hamilto- 
nian not supporting any coherent generation of net RP 
nuclear polarization, as in the example considered pre- 
viously, it will be (I z ) = Tr{pl z } = 0. This necessarily 
implies that {I z ) s — — {I Z ) T - To parenthetically prove 
this, note that since I z commutes with the electron spins 
Si and S2, it also commutes with Qs = 1/4 — si ■ S2 and 
with Qt = 1 — Qs- Thus, since the two projection oper- 
ators are orthogonal, i.e. QsQt = QtQs — 0, it will be 
QsIzQt = QtIzQs = 0. Furthermore, since Qs + Qt = 
1, it will be Tr{pL} = Tr{p(Q s + Q T )I z (Qs + Qt)}, 
and by using the previous relations it follows that indeed 
= (I,) = (I,) s + (I,) T . 

Now, the main idea behind the new kind of CIDNP 
signal generation is this: the intramolecule S-T deco- 
herence is equivalent to some RPs being projected to ei- 
ther the singlet or the triplet state, exactly due to the 
continuous quantum measurement of Qs induced by the 
RP recombination dynamics, the measurement rate being 
(ks + kx) /2 [23|. The result of this measurement is either 
Qs = 1, i.e. the RP is projected to the singlet manifold, 
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or Qs = 0, i.e. the RP is projected to the triplet mani- 
fold. The probability for a singlet (triplet) projection is 
(Qs) ((Qt))- Hence, at time t+dt the state of the radical- 
pairs will be (ignoring the RPs that recombine into neu- 
tral products during dt) a new mixture comprised of (a) 
radical pairs in the state p + dp, with dp = —idt[H,p], 
plus (b) the singlet-projected radical pairs described by 
QspQs, plus (c) triplet-projected radical pairs described 
by QtpQt- The nuclear polarization of states (a) will 
still be zero, since again, the Hamiltonian % does not 
generate net nuclear polarization. The combined nuclear 
polarization of projected states (b) and (c), denoted by 
(I z ) pT °i, will be the weighted sum of singlet RP nuclear 
polarization, (I z ) s , and triplet RP nuclear polarization, 
(I Z ) T , the sum being weighted by the respective projec- 
tion probabilities, i.e. 
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Taking into account that (I Z ) T = — (I z ) and using the 
completeness relation (Qt) = 1 — (Qs), it follows that 



(I z r°i = (I z ) s (2(Q s ) 



- 1 



(2) 



Thus, if (Q s ) + (Qt) + 1/2, more (when (Q s ) > 1/2) 
or less (when (Qs) < 1/2) projections will take place to 
the singlet state than to the triplet, and (7 z )f IOJ will be 
non-zero, if of course (I z ) s is non-zero. 

We can now estimate the quantum measurement cor- 
rections to CIDNP, to be denoted by (I z ) qc , shown in 
Fig. 2D. Letting k = k$ = kT, the rate of singlet or 
triplet projections (the "quantum measurement" rate) 
is just (ks + fcr)/2 = k. We will focus on low mag- 
netic fields and assume, as is the case most often, that 
k > A, lu. Towards the estimate, consider the following 
three steps: (1) As the RPs start out in the singlet state, 
and the S-T mixing frequency ft <C k, in the short re- 
action time of 1/k <C 1/ft the singlet character will not 
have enough time to change appreciably, so it will be 
(Qs) ~ 1- (2) The electron spin of the donor molecule, 
which contains the one and only nuclear spin of the RP 
under consideration, will precess at frequency \lj ± A\ 
when the nuclear spin is in the |±) state. So at time t, 
the number of singlet RPs with the nuclear spin in the 
±) state will be N± « cos[(w ± A)i] 2 . The nuclear po- 
larization of singlet RPs at time t = 1/k will thus be 
(I z ) s = (N + - N-)/(N + + N-), easily seen to be 
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(3) 



And (3), S-T projections taking place at the rate k will 
be a source term for d(I z ) qc /dt. The source term, pro- 
portional to (7 z ) pr °j, starts at t = from zero and 
oscillates at a rate ft while it is drained at a rate k, 
thus d(I z ) qc /dt = fc(/ 2 )P roj e- fc *sinf2i, yielding at time 
t w 1/k that (I z ) qc w -(I z )v ro m 2 /k 2 . Putting every- 
thing together and using @, we finally arrive at 
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For the particular example depicted in Fig. 2 it is (Fig. 
2B) ft « A/10, hence according to the estimate (|4j) it 
is (I z ) qc ~ 4 x 10~ 6 , in good agreement with the ex- 
act value shown in the right y-axis of Fig. 2D. Further- 
more, the result (g]) can be recast in terms of the thermal 
polarization 7 t h = fuxj/A-yksT (the expression for / t h is 
valid for ftw <C ^ksT), where w = "f e B is the electron 
Larmor frequency in the external magnetic field B, with 
7e = 27rx 2.8MHz/G being the electron gyromagnetic 
ratio and 7 = \p e /p p \ = 658.5 is the ratio of electron- 
to-proton magnetic moment. The enhancement factor 
(I z ) q c/ h\\ is thus seen to be independent of uj and equal 
to 
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(5) 



For ease of use, in the previous expression we provided 
the enhancement factor in terms of the S-T mixing fre- 
quency fl (in units of 0.01 ns -1 ), the hyperfine coupling 
A (in units of 0.1 ns -1 ) and the recombination rate k 
(in units of 1 ns -1 ). For typical hyperfine couplings of 
about 10 Gauss and recombination times on the order 
of 1 ns, it is seen that the quantum measurement correc- 
tions to chemically induced dynamic nuclear polarization 
amount to a significant enhancement of at least four or- 
ders of magnitude over the thermal polarization. 

We note for the reader who would worry about angular 
momentum conservation that the RP's nuclear polariza- 
tion is exactly balanced by an opposite and equal elec- 
tronic polarization, i.e. the sum (si z + S2 Z + I 2 ) = at 
all times, as it should be, since it is zero at time t = 
(the initial state is po — Qs/Tr{Qg}, which is a sin- 
glet state with zero nuclear polarization). The non-zero 
electronic polarization necessarily involves the admixture 
of the \T + ) and |T_) triplet states, therefore this new 
CIDNP mechanism is a low-magnetic-field mechanism, 
because at high fields the \T±) states are out of reach 
from the singlet state \S). How low a magnetic field de- 
pends on the recombination rate fc, since the triplet RP 
decay broadens the triplet energy levels by about k, hence 
the mechanism is appreciable for magnetic fields B such 
that B < k/j e , i.e. when the Zeeman energy separation 
of \T±) is within the broadening k. For recombination 
rates on the order of 1 ns -1 , the effect is appreciable for 
fields B < 50 G. 

In closing, the following comments should be made: 

(1) Kaptein's rules provide a compact relation of the 
main features of the CIDNP signal to the system's pa- 
rameters. For this case, the sign of the effect, i.e. whether 
it is emissive (+) or absorptive (-) is the same as the sign 
of the hyperfine coupling A. 

(2) The new CIDNP effect reported here does not re- 
quire any particular combination or fine-tuning of mag- 
netic interactionSjhke the coherent TSM low-field ef- 
fect analyzed in [2J], but is rather general and shows 
up whenever there is a non-zero nuclear spin polariza- 
tion of singlet RPs, (I z ) s , as in Fig. 2B. The rather 
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general requirement for this to happen is that the oscilla- 
tion frequencies uj m — uj n of the non-zero matrix elements 
[Qs]mn = (w|Qs|n) must have terms linear in the hy- 
pcrfine couplings. Here |n) and uj n are the eigenvectors 
and eigenvalues of H, with n = 1, 4M, where M is the 
nuclear spin multiplicity of the radical-pair. 

(3) Experimentally, while NMR measurements at low 
fields suffer loss of sensitivity, recent methods of ultra- 
sensitive magnetometry using, SQUIDs dHHH or atomic 
magnetometers [HHHf, have demonstrated the ability to 
sensitively detect NMR signals at zero magnetic fields. 
So the detection of the new CIDNP signal analyzed here 
is rather straightforward. 

(4) Although we considered the special case ks = fcr, 
the new kind of CIDNP signal we predict will ubiqui- 
tously affect all sorts of radical pairs, including of course 
the ones with asymmetric recombination rates ks ^ kx, 
which predominantly occur in photosynthetic reaction 
centers [2l[ . We have chosen to deal with the special case 
ks = kx and a magnetic Hamiltonian unable to support 
coherent generation of CIDNP signals in order to be able 
to focus exclusively on this fundamentally new CIDNP 
process. Establishing how this mechanism attributed to 
S-T decoherence affects, qualitatively and quantitatively, 
the cases where there exists a CIDNP signal according 
to the mechanism's traditional understanding (e.g TSM, 
DD or DR scenarios) is more involved an exercise that 



will be undertaken elsewhere. 

(5) We have considered the CIDNP signal in the RP 
state, and not in the diamagnetic ground state of the 
molecule in a closed reaction cycle as in Fig. IB. How- 
ever, it is easily seen that the same results hold for the 
ground state in the steady state under continuous exci- 
tation and a closed reaction cycle. 

Concluding, this is the first step of a very promising 
research program, namely what are the particular reper- 
cussions of the recently unraveled fundamental quantum 
dynamics of the radical pair mechanism for understand- 
ing the dynamics of photosynthetic reaction centers. We 
have extended the fundamental understanding of CIDNP 
by pointing to a new and qualitatively different way to 
obtaining enhanced nuclear spin polarization. Interest- 
ingly, this polarization naturally turns out to be signif- 
icant at low magnetic fields pertinent to natural photo- 
synthesis. It stems from the random singlet or triplet 
projections taking place during the quantum evolution 
of radical pairs. S-T projections are the source of S-T 
decoherence fundamentally present in RP recombination 
reactions. Thus, an intriguing question raised is if and 
how the quantum mechanical process of spin decoherence 
affects the biological workings of photosynthetic reaction 
centers. In other words, what is the biological significance 
of earth-field CIDNP produced by S-T decoherence? 
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